Slow Light Propagation in a Thin Optical Fiber via Electromagnetically Induced 

Transparency 
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We propose a novel configuration that utilizes electromagnetically induced transparency (EIT) 
to tailor a fiber mode propagating inside a thin optical fiber and coherently control its dispersion 
properties to drastically reduce the group velocity of the fiber mode. The key to this proposal is: 
the evanescent- like field of the thin fiber strongly couples with the surrounding active medium, so 
that the EIT condition is met by the medium. We show how the properties of the fiber mode are 
modified due to the EIT medium, both numerically and analytically. We demonstrate that the 
group velocity of the new modified fiber mode can be drastically reduced («44 m/sec) using the 
coherently prepared orthohydrogen doped in a matrix of parahydrogen crystal as the EIT medium. 

PACS numbers: 42.81. Dp, 42.50.-p, 42.50.Gy 



I. INTRODUCTION 

There has been many excitements on the group veloc- 
ity management via electromagnetically induced trans- 
parency to drastically reduce the group velocity of light 
[0, ^ ||, to freeze to store the quantum infor- 
mation of light to achieve superluminal group ve- 
locity propagation and to control the group velocity 
between subluminal to superluminal propagation |^ in 
varieties of resonant mediums. Many experiments have 
been reported on reduction of group velocity of light 
(slow light) using ultra-cold atoms |l], in hot atoms 0] 
and also recently in a rare-earth ion doped crystal ^f. 
Slow light has also been observed in far-off-resonance 
Raman systems Extensive study of the slow light 
can be found in many recent literatures Ml including 
many interesting applications ||l^, |ll|, |l^, . Harris 
and coworker have demonstrated possibility of observing 
non-linear processes at low light level using the slow light 
[0. Matsko et al. have demonstrated enhancement of 
acousto-optical effects in a dielectric fiber doped with res- 
onant three-level system when velocity of light becomes 
comparable to velocity of sound ||T^. Possible applica- 
tions of slow light for efficient four- wave mixing 131 and 
polarization splitting Q are also reported. It has also 
been observed that dynamically turning "ofF and then 
turning "on" the strong control field, a signal pulse can 
also be freezed and retrieved ||, ^, leading to possible 
realization of a quantum storage device. The key to all 
the above phenomena is existence of large dispersion be- 
havior combined with the very low absorption due to EIT 
[Q . The most attractive feature of the EIT mechanism 
has been the dynamical control of the dispersive property 
of the medium using the control laser parameters. 

In recent years rapid progress has also been made in 
slow light propagation in photonic heterostructures. It 
is reported that extremely slow light can be observed in 
photonic band gap materials ||l5|, in Moire fiber Bragg 
gratings [ p^ , and in optical resonator-array waveguides 
[ jl7[ . In these kind of studies the material dispersion, 



which is determined by the engineered structure, is used 
to achieve the slow light. These materials have significant 
importance due to their applications to scalable photonic 
devices. However the structure dependent dispersion- 
management lack in the dynamical control compared to 
the EIT assisted control of dispersion. 

In a conventional optical fiber, a propagating single- 
mode field experiences a dispersion due to the dispersive 
material contents of the core and clad. The group veloc- 
ity of the envelope of such a propagating fiber mode is 
given by 
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(1) 



where [3 is the propagation constant of the fiber mode 
and u is the angular frequency. In conventional opti- 
cal fibers such dispersion is very small compared to the 
dispersions in fibers using tailored materials such as pho- 
tonic heterostructures |15, ^ pT[ . 

In this paper we propose a configuration to tailor the 
dispersion properties of a fiber mode propagating in- 
side a thin (tapered) optical fiber surrounded by an EIT 
medium [see Fig. ||], and as a result we show that the 
group velocity of the fiber mode can be reduced signifi- 
cantly. Our method utilizes EIT and hence has the ad- 
vantage of large control over the dispersion property of 
the fiber mode. It should be noted that our usage of 
"thin fiber" should be understood as the case when the 
diameter of the fiber 2a ~ A, A is the central wavelength 
of the field that propagates inside the fiber. Under such 
a condition large part of the energy (evanescent-like) of 
the fiber mode lies outside the thin fiber. Thin tapered 
fiber has been a standard method in the dielectric sphere 
experiments to get large coupling of the output field with 
the whispering-gallery modes of the sphere p8[ . 

Further we have introduced ortho-H2 doped in a ma- 
trix of para-H2 as an EIT medium surrounding the fiber 
to demonstrate the above mentioned reduction in Vg. 
This medium has the following advantages over other 
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EIT mediums: (a) It has very high density and short 
interaction length (few tens of ^m) compared to its gas 
phase, (b) The vibrational states considered for the EIT 
transitions have very small dephasing (^53 Hz at an 
operating temperature of 4K To our knowledge, 

the only other solid in which EIT has been reported is 
Pr:YSO crystal Using our scheme we achieve a 

large dispersion of (3 and hence a huge reduction in the 
group velocity (~44 m/sec for a chosen set of field pa- 
rameters). We also predict the possibility of stopping the 
fiber mode. With the rich technologies available for the 
optical fibers, we hope, slow light in optical fiber could 
open up new possibilities in non-linear optics studies and 
for integrated photonic device applications. 

The organization of the paper is the following: In Sec. 
II, we present the propagation of a fiber mode inside 
a thin single-mode optical fiber, particularly when the 
medium surrounding the thin fiber is an active medium. 
In Sec. Ill, we describe how the behavior of the fiber mode 
can be modified when the surrounding medium is an EIT 
medium and introduce a dressed fiber mode. In Sec. IV, 
we obtain an approximate analytical expression for the 
group velocity of the new modified fiber mode and show 
how the group velocity can be substantially reduced by 
choosing the appropriate parameters. In Sec. V, we give 
a detailed account of realization of slow group velocity of 
the fiber mode using ortho-Il2 doped parahydrogen crys- 
tal as an EIT medium surrounding the fiber. We present 
more detailed aspects of the slow light propagation of the 
dressed fiber mode with numerical results in Sec. VI. We 
summarize our results in Sec. VII. 



II. PROPAGATION DYNAMICS INSIDE A 
THIN FIBER SURROUNDED BY AN ACTIVE 
MEDIUM 

A thin fiber is obtained by heating a conventional com- 
munication fiber and pulled to obtain a narrow waist . 
The waist of thin fiber is 2a ~ A. Thus it has almost a 
vanishing core and hence the refractive indices that de- 
termine the fiber modes are the cladding refractive in- 
dex which we denote as n/ (and will be referred as fiber 
refractive index henceforth) and the medium refractive 
index rim surrounding the thin fiber. 

The general equation for propagation of an electric field 
Ect through a medium with a refractive index n is gov- 
erned by 
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(2) 



where, c is the velocity of light in free space, and a is a 
field index that distinguishes the fields, when more than 
one laser field propagate through the fiber. Let us con- 
sider propagation of a pulse inside the thin single-mode 
fiber. A pulse corresponding to the lowest order charac- 
teristic mode of the fiber with central frequency uja may 
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FIG. 1: The configuration under consideration: a thin 
tapered optical fiber is surrounded by the EIT atoms or 
molecules. Two fields, a strong and a weak field, prop- 
agate through the single-mode thin optical fiber. The 
evanescent-like part of the fiber modes interact with the 
EIT atoms/molecules within the interaction region as shown 
above. The energy- level scheme for the EIT atoms/molecules 
is shown in the inset. The spontaneous decay coefficients are 
denoted as 27^, the strong control field and weak probe field 
Rabi frequencies are given by 2G and 2g respectively, their 
corresponding frequency detunings, from the transitions they 
couple, are A and 5. The ground state dephasing is denoted 
as r. More details of the EIT scheme is described in Sec. III. 



be given by 

£:„(r,0;t) =4 / £a(r;w)e*''°^-*"*dw + c.c. (3) 



Here ^^(r;^) is the distribution of the field amplitude 
at frequency lu in the transverse direction and Pa corre- 
sponds to the propagation constant of the propagating 
mode with uj^ being the central frequency. On substi- 
tuting (H) into (H) the equation for the field amplitude is 
obtained as 



9y2 



(4) 



where ka = Wq/c is the free space propagation constant 
corresponding to Ua. The refractive index n = Uf (rim) 
inside (outside) the fiber. In deriving Eq. (||), we have 
made the following approximations: (a) The spatial vari- 
ation of the refractive index along the transverse direc- 
tion is assumed to be zero inside the fiber and is negligible 
outside the fiber (i.e. Vn ~ 0). (b) The time-duration of 
the pulse is assumed to be long compared to the response 
of the medium, so that the medium can be considered to 
be in steady state {drim/dt = 0). (c) Further we work 
in the single (lowest order) mode of the fiber - i.e. under 
the condition [E2| 



A. 



27ra 



< Ao 



(5) 
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Here Aq = 27r/fcct, a is radius of the thin fiber, (c = 
1.405 corresponds to the first zero of the Bessel func- 
tion Jo(Cc) = 0. (d) We also use the weakly guiding ap- 
proximation {{uf — nm)l_nm ^ 1) for which the approx- 
imate solution of Eq. (0) corresponding to the lowest- 
order mode HEii is a linearly polarized LPqi mode 

The propagation equation (jij) for fiber has well-known 
analytical solution for the fundamental mode Pq, p2| 



""°/"/"^°^\ for.<a 



for r > 



(6a) 



(6b) 



which is obtained by matching the boundary conditions 
for normal and tangential components of the field sep- 
arately and is normalized. Here r = \/ x"^ + is the 
radial distance from the axis of the fiber, a is the radius 
of the thin fiber and is a constant that determines 
the amplitude of the field, and 
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(7) 



The symbols Ji and Ki represent Bessel functions of first 
kind and modified Bessel functions respectively of order 
I. In writing the above solution the field outside the 
fiber is approximated to an exponential for case of 
further evaluation of the integrals. The characteristic 
equation that determines the propagation constant (3a 
for the fundamental mode is 
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(8) 



The radial distribution of the energy of a fiber mode 
propagating inside the thin fiber, given by Eq. (|^), 
clearly depends on the frequency of the propagating field, 
fiber parameters such as diameter of the fiber (2a), re- 
fractive index Uf, and refractive index of the medium Um- 
For example, for kaa = 1.31 and nf = 1.43, about 49% 
of the total energy of the fiber mode lies out side the thin 
fiber when refractive index of the medium is — 1. 

Generally the medium outside the fiber is air or a 
passive medium having constant refractive index n™, 
for which the above mode calculation is valid. But in 
the present case under consideration, we have an active 
medium that exists surrounding the fiber, where the re- 
fractive index of the medium n„i = 7i,„(r) is no more 
a constant, moreover it could be a function of parame- 
ters of the propagating field. In such a case the above 
description of mode solution fails. However we use an in- 
direct and approximate approach to determine (3a using 
an effective refractive index for the medium to replace 
rim, which is described as follows: 



Let us introduce an effective refractive index Um for the 
active medium outside the thin fiber, experienced by the 
field £p. The corresponding polarization of the medium 
at a radial distance r is given by 



Pa^^ir;uja) = eo('^m " l)£a{r;UJa)- 



(9) 



Here eo is the electrical permittivity of the vacuum. How- 
ever the actual polarization of the medium has a radial 
distribution, given by 

Pa,^{r;uJa) = eQ{nl^{r) - l)£air;uja)- (10) 

For rim and hence to represent the actual physical 
system described by Eq. (|l0|), we demand that the total 
interaction energy densities £* ■ Pa^,^ and £* • Pa^ over 
the whole transverse plane due to Pa^v and Pa^ , respec- 
tively, should be the same. Thus multiplying £*{r]uja) 
to both Eqs. (^ and (|lO|), integrating over the transverse 
direction, and equating their right hand sides we get 



la £a{'r)nl,{r)£a{r)rdr 
l^\£airWrdr 



(11) 



in cylindrically symmetric coordinates. Here R is the ra- 
dial extent of the medium from the axis of the fiber. Dif- 
ferent schemes for averaging refractive index are used to 
determine the field propagation in heterogeneous struc- 
tures (e.g. see in pi]). 

Therefore n„j in Eqs. (^||) can be replaced by the aver- 
aged refractive index of the medium hm that is constant 
along the radial direction. Thus parameters in Eq. (|^) 
are modified to 



«/c = ^Jklnj- (31 

ipa = kaCL^Jri' 



2 

m t 



, (12) 



that determine the field distribution £a (r) and the prop- 
agation constant (3a of the fundamental mode. However 
it should be noted that nm{r) in Eq. ( pT|) itself is a func- 
tion of ^^(r). Therefore, the propagation is initiated, for 
example, with a Gaussian field at the input, and soon 
that evolves to the fundamental mode of the fiber £a{r) 
given m Eqs. (1) and (|l|). 



III. DRESSED FIBER MODE VIA EIT 

In the previous section we have shown that the 
evanescent-like part of the fiber mode can be coupled 
with the medium outside the fiber for suitable sets of 
field, fiber and medium parameters. And thus if the 
medium outside the thin fiber is active, the fiber mode 
will closely follow the nature of the active medium. In 
this section, we demonstrate how a fiber mode can be tai- 
lored using an EIT medium surrounding the thin fiber. 



4 



Let us consider that the active medium consists of 
low density homogeneously broadened atoms / molecules 
having three-level A schemes for their energy levels (see 
the inset of Fig. |l]). We assume the atoms / molecules 
do not move and hence we ignore the transit time broad- 
ening. The |1) ^ 1 2) and |1) ^ |3) are coupled by a 
strong control field Ec{r) and a weak probe field Ep{r) 
respectively. The control field fiber mode (CFM) and 
probe field fiber mode (PPM) are given by Eq. (^) with 
a — > c for the control field and a p for the probe 
field. The |2) ^ |3) is a dipole forbidden transition. 
We assume that both PPM and CFM satisfy the single- 
mode condition (|^). The refractive index of a general- 
ized EIT medium, corresponding to a weak probe field 
is well known Therefore in the following, we only 

present the formula for refractive index in our configura- 
tion, where the constant control field Rabi frequency in 
10 is replaced by the transversely distributed function 
2G(r) = 2d£c{r)/fi; where d is the induced dipole mo- 
ment corresponding to |1) <-!■ |2) transition, h is Planck 
constant. Refractive index of such a medium experienced 
by the weak PPM £p in presence of the strong CPM £c{r) 
can be obtained as 



i-fi[r - i{A - S)] 



ill + 72 + «'5)(r 



z(A-(5))-K|G(r)|2 

^ (13)' 
for small absorption. In deriving Eq. (O) we have used 
following assumptions: (a) the two field modes are un- 
coupled and do not interact with each other during the 
propagation inside the fiber, (b) the CPM £c propagates 
without loss inside the fiber. We will present a more de- 
tailed derivation of the refractive index, for our configura- 
tion, using a specific example in Sec. V. It may be noted 
that the medium contribution to the refractive index for 
the control field is negligibly small as the control field acts 
on |1) 1 2) transition where very small population or no 
population is available to make any significant contribu- 
tion to the susceptibility, and hence the refractive index. 
Thus for CPM propagation rim = — 1, using vacuum 
as the medium outside the thin fiber. This n,„ is used to 
calculate the fiber mode for the control field £c- Here the 
detuning of the central frequency of PPM (CPM) from 
the transition frequency ujq of |1) ^ |2) {ujq of |1) <-> |3)) 
is given hy S = ujQ — ujp (A = — Wc); '^li represents the 
spontaneous decay coefficient corresponding to |1) \i) 
transition and 2r represents the dephasing rate of the 
ground state. It may be noted that for our configuration 
in Fig. 0, the control field Rabi frequency 2G = ^.Gir) 
and hence the refractive index Um = nm{r) are functions 
of r which is different from standard EIT configurations 
[Q. In Eq. (|l^) we have used a characteristic (dimen- 
sionless) parameter ^ for the EIT medium that is defined 



as 



rieo7i 



(14) 



Here N is the density of atoms / molecules. For typical 
EIT experiments with alkali atomic gases (see references 



in jlj]) ^ = 0.01 ^ 0.1. More detailed calculation on a 
specific realistic system is presented in Sec. V. 

Assuming the susceptibility of the medium Xm to be 
small and denoting h as the average refractive index over 
the whole transverse cross section of the fiber, the prop- 
agation constant /3p for the fiber mode of the probe field 
can be approximated as 



C J^^ £*n£prdr 
[ C\£p\'rdr , 
{brim + (1 - b)nf} kp. 



(15) 



Here rim, defined in Eq. 
ten as 



\/nj , can be approximately writ- 



i:\£p\'rdr 



(16) 



and the coefficient b denotes the fraction of total energy 
of the fiber mode outside the fiber and is given as 



lo \£p\^rdr 



(17) 



Clearly from Eq. (|T|), Pp of the PPM will be modified by 
the medium for any non-zero value of b. Next we proceed 
to calculate b. 

The distribution for £p{r), for r < a, is governed by 
Bessel function Jo{i^fpf) as in Eq. (|6^). We expand 
~ V4 , (4p^V4)^ 



■ — . . . and consider only 



^ (2!)2 

the first two terms in the expansion for the thin fiber 
where Kf^r is small. Thus substituting £p from Eq. (|6 
in (frh and carrying out the integration we obtain 



1 



Jp 



08)" 

where <j)p and k/ are the coefficients defined in Eq. (|l^). 
Here we consider the extent of the medium in the radial 
direction is very large, i.e. R —> oo. It can be easily 
shown that 5=1 for a — > (and hence (j)pa, Hf^a — s- 0). 
Therefore the smaller is the radius a, larger is the effect 
of the surrounding medium on j3p. For the parameters 
used in Fig. H, numerical value of 6 0.57. 

In Fig. H we present the dispersion and absorption 
characteristics of the new modified PPM. In Pig. ^a), 
we plot the solution (scaled with fco) of the charac- 
teristic equation Eq. (H). The numerical method for the 
plot is described in Sec. VI. Clearly the PPM assumes the 
same structure of dispersion and absorption as that of the 
EIT medium and (3p = l3p{G, A, S). The dashed lines in 
Fig. 1^ represent dispersion and absorption characteristic 
of the weak PPM when CFM is absent. In presence of 
the CFM, the absorption of PPM becomes very small at 
S = 0. We name this EIT modified PPM with propa- 
gation constant Pp as dressed fiber mode. The origin of 
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FIG. 2: (a) The dispersion and (b) the absorption characteris- 
tics of the PFM. The dashed hnes correspond to the behavior 
of the PFM alone, and the soHd Unes depict the modification 
occured in presence of the CFM (i.e. due to EIT). Clearly 
the behavior of the PFM is governed by the properties of 
the EIT medium. The parameters used here are: ^ = 0.107, 
71 = 72 = 7 = IMHz, r = 0, n/ = 1.43, 2a = 0.3/im, 
A = 780nm, £c{r = 0) = h"f/d and A = 0. All frequencies are 
scaled with 7. 

this new behavior of the fiber mode can be understood as 
due to the strong coupling of the evanescent-like part of 



the fiber modes with the surrounding EIT medium. This 
is the key for the next sections and for this paper which 
helps us realize slow light in the thin optical fiber. 



IV. SLOW GROUP VELOCITY OF THE FIBER 
MODE 

In this section we examine the group velocity of the 
propagating new dressed fiber mode described in the pre- 
vious section. For PFM with central frequency tOp — loq 
which is on resonance with the |1) <-> |3) transition and 
propagating in the single-mode optical fiber, the group 
velocity in Eq. p) can be rewritten as 



d(3p 


-1 


dujp 





It is now clear from Fig. [^(a) that the steep dispersion 
of (3p around cjq can cause a large reduction in group 
velocity. 

To obtain a clear physical picture of the contributions 
from various dispersions to the group velocity of the fiber 
mode, we use the approximate relation ( |l5| ) and substi- 
tuting into Eq. ( p^ ) we get 

^-^^(^^™ + (l-%/)U.., (20) 

Dropping the dispersion associated with the fiber-clad 
[dn f /dojp] that is negligible compared to that of the res- 
onant medium [dnm/dcup] and substituting from Eq. 
dll) we get 



UJQ^ 

c 



, db b 



dojp ^ 



Elrdr + 2 



(21) 



In deriving the above equation £p is taken to be real. We 
have also dropped the terms associated with dkp/dojp 
that are negligibly small (~ 1/c) . From Eq. ( p^ we get 



dun 



duj„ 



71 f 



|G(r)| = 



[\G{rW + hi+l2)rY 



(22) 



at the EIT condition and when the control field CFM 
is on resonance with |1) |2) transition. Further us- 
ing (6b) for £p{r) and Sdr) and using (p2|), we obtain 
an expression for the group velocity of the PFM in the 



presence of the CFM 

1 cjo7i? ' 



{1 + 2(0p - (j),)a} 



2cGl 

Wo 



(j)cY{l + 2(l)pa) 
^ db 

duJri 



(23) 



where Go = G{r = a). We have considered that the 
medium is present over large radial distances (i? — > 00). 
The parameters (pp and are the PFM and CFM param- 
eters defined in Eq. ( [l2|) with a ^ p and a ^ c respec- 
tively. We note that in deriving (^3|) we have neglected 
the third term in ( [2l| ) compared to the other two terms, 
as it contains the small multiplicative term («„ — n^). 
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Numerically we observed that the third term is at least 4 
to 5 orders of magnitude smaller compared to the other 
two in (^) , for a thin fiber surrounded by a typical EIT 
medium. We have derived the above expression for group 
velocity at the ideal EIT condition, i.e. for F = 0. It is 
interesting to note from Eq. ( p3|) that, assuming the sec- 
ond term is small compared to the first, the fiber mode 
can be brought to a complete halt {vg = 0) by turning off 
the CFM (Go = 0), which could be useful for storage |^ 
of light information in the PFM. Such dynamical aspects 
of the group velocity of the fiber mode will be discussed 
elsewhere. 

Further we examine the limiting case of vanishing fiber 
radius (a ^ 0), where (as 13c ~ kc, fim — rim — 

1) and hence 0c ^ 0. Also 6 ^ 1 as discussed in the 
previous section and hence db/dup ~ 0. Substituting 
these in Eq. (|2^ ) we get 



(24) 



Thus we recover the expression for group velocity of 
a pulse propagating inside a bulk sample of the EIT 
medium with the zero ground state dephasing, in the 
case of vanishing fiber radius. 



V. ORTHO-iJa DOPED IN SOLID PARA-Ha 
CRYSTAL AS AN EIT MEDIUM 

We propose the orthohydrogen molecules (with nuclear 
spin 1=1) doped in the solid parahydrogen (/ = 0) 
molecular crystal as an EIT medium. The transition un- 
der consideration is a crystal field induced vibrational- 
rotational transition Qi{l) (w = 1 <— 0, J = 1 ^ 1, AIj = 
— 1 ^ 0, ^ 1). The degenerate sublevels of rotational 
angular momentum J = 1 states with mj = 0,±1 split 
in the anisotropic crystal field. Further the quadrapole 
moment of the orthohydrogen (J = 1) molecule in- 
duce dipolc moments on the surrounding parahydrogen 
(J = 0) molecules. And these many body dipolc mo- 
ments interact with a radiation, having its electric field 
perpendicular to c-axis of the crystal, to cause the transi- 
tion V = V 



1 



26, 2711. This transition occurs at 



infrared Aq = 2.4/im and the selection rule for this tran- 
sition is Am J = ±2 [See Fig. |^(a)]. Further applying 
an axial magnetic field, the degeneracy in the hyperfine 
levels 1=1 can be removed to obtain a six level scheme 
as depicted in the Fig. ^(b). 

Before we proceed, we note that the selection of this 
medium as the EIT medium is of particular importance 
to the present problem, due to the following: (a) As dis- 
cussed earlier we need the field to satisfy single mode 
condition (||) and also a large evanescent-like part of 
PFM should be available for interaction of the fiber mode 
with the medium, for which 2a < Aq. Thus for ortho-H2 
medium (Aq = 2.4/Ltm) one can work with 2a ~ 1 — 2fim - 
whereas for a standard EIT medium having optical tran- 
sitions, the thickness of the fiber has to be as small as 




FIG. 3: (a) The crystal field induced = <-> w = 1 transi- 
tion. The selection rule for this transition is Amj = ±2. (b) 
The six level scheme that can be obtained from the hyperfine 
splitting of 7 = 1 degenerate levels applying an axial mag- 
netic field. Here 27^ are the non-radiative decay rates from 
1; = 1 — > n = states, F is the crystal field induced mixing 
rate of the ground states, 2Gi and 2gi are CFM and PFM 
Rabi frequencies defined in Eqs. (^i]) and (^), respectively, 
(c) The effective three-level A scheme reduced from the 6-level 
scheme which describes the EIT in the ortho-Ha molecule. 



2a ^ 1/im. (b) This solid has very narrow widths |2^ 
compared to the other solid EIT medium j|] . (c) Further 
we note that since we are working here with the solid 
medium, the assumption of neglecting the transit time 
broadening in Sec. Ill is easily met. 

The magnetic degenerate sub-levels in Fig. ||(a) can be 
coupled selectively by choosing suitable polarizations of 
the laser fields Ig^l- We use a strong a+ polarized laser 
field that couples |2) ^ |4) and |3) <-> |5) transitions, 
and a cr_ polarized weak probe field that couples the 
|1) ^ |5) and |2) |6) transitions. Thus the interaction 
Hamiltonian in Rotating Wave Approximation (RWA) 
can be written as 

-Hi = -;i[(.gi|2)(6|+g2|l>(5|)e-*'^''*+'^-^ 

+ (Gi|3)(5| + G2|2)(4|)e-'"=*+"3=" + H.c] (25) 

Here ujp (cUc) represents the central frequency of the PFM 
(CFM), and (3p (/3c) represents the propagation constant 
of the corresponding propagating mode. The Rabi fre- 
quencies of the CFM are. 



(jl — r , (-r2 — — 



and that of the PFM are 



91 



92 



h ' 



\di5\£p 



(26) 



(27) 



Here dij = {i\d\j) represents an effective dipole moment 
corresponding to the crystal field induced |z) <-!■ \j) tran- 
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sition, which is calculated from the experimentally mea- 
sured value of the integrated absorption coefhcient [p5| . 
The unperturbed Hamiltonian in terms of Fig. ^(b) can 
be written as 

Ho = + 217)|1)(1| + +17)|2)(2| 

+?iw„|3)(3| + 2ar2|4)(4| +ni7|5)(5|. (28) 

Here Uoov is the separation between the states |n) and \n-\- 
3) (for n = 1,2,3) in Fig. |(b). The hyperfine splitting 
of / = 1 degenerate levels, caused by the axial magnetic 
field, is represented by 251. 

Thus the equation of motion for the six-level system is 
given by the density matrix equation 

3 r, 3 6 



calculations presented below we have taken all the six 
levels into consideration. The conservation law for the 
population in different energy levels is 



(32) 



7. { N) I , P}+ + 3 I J> (j I 

1=1 i=l j=4 



i=4 



(29) 



Here 27iS could be the non-radiative decay rates from 
each of the v = 1 ^ v = states and 2r is the crys- 
tal field induced mixing rate of the ground states. The P26('") = 
notation { }_|_ represent the anti-commutators. Thus the 
equations for the density matrix elements are 

P22 = -272P22 + iG2P42 - iG*2P2i + i.g2P62 " iglhfi 
P2i = -(72 + 2r + z(A - ri))(524 + «G2(P44 - P22) + i.giP64 

P26 = -(72 + 2r + i((5 + f7))p26 + iG2P46 + i5i(p66 - P22) 

P44 = -4rp44 + 2r(p55 -I- Pee) 
2 

+ 3 (71PII + 72/022 + 73P33) + iG*2P2A - iG2PA2 

P46 = -(4r + i(5 - A + 20))p46 + iG2P26 - «5lP42 
P66 = -4rp66 + 2r(p44 + P55) 



By the particular choice of the field configuration (see 
Fig. |4(b)), the strong CFM Gi{r) pumps the population 
from^4) and |5), and via the non-radiative decays from 
1 2) and |3), a good initial state |6) can be prepared in the 
molecules close to the fiber. From (^) we numerically 
estimated the steady state value of ground state pee {i") ~ 
0.97 when the CFM Rabi frequency at r is G(r) = 7, A = 
0. Here we have assumed that 7^ = 7 (i = 1, 2, 3) and we 
have used the parameters 27 = 30 kHz |2^ and 2r = 53 
Hz The Rabi frequencies are assumed as \Gi\ = 

G (i — 1,2) for simplicity. Thus the configuration in 
Fig. U(b) effectively reduces to a three-level A scheme as 
shown in Fig. ^(c) in presence of the cr-(- polarized strong 
CFM (Gi ^ 0), and hence the dynamics is primarily 
governed by the Eqs. (^). 

Assuming that the steady state population of |6) in 
absence of PFM is equal to unity, we get the solution of 
the density matrix equation for a weak PFM as 



igi{r)[4:r + i{S - A + 2n)] 



(7 -f 2r -f i{s + n)){4r + i{s~A + 20)) + |G(r)|2 • 

(33) 

The frequency dependent susceptibility of the medium 
containing orthohydrogen doped in parahydrogen crystal 
can be written as 



Xpa 



^<^26, 



(34) 



where Xpaia is susceptibility of the background parahy- 
drogen matrix, a26 = 7P26/gi and ^ = N d^g / (Heoj) 
(where d^s — dij effective dipole moment, assuming the 
effective dipole strengths of all involved transitions to be 
same). All other parameters are same as defined earlier. 
Thus the refractive index of this doped medium is 



+ 2 (71P11 + 72/022 + 73/O33) + iglhe - mP62- 



(30) 



n,n{uj) = (1 + Re x™)^/2 



e(Re (726)]'/', (35) 



In writing Eqs. (|3C|) we have eliminated the temporal 
and spatial rapid oscillations by making the following 
transformations 

P24 ^ p24:GXp[-iLUct + iPcZ], 

P26 P26exp[-zcjpi 4-i/3pz], 

P46 P46exp[-z(cJp-tJc)i + «(/3p-/5c)2], 

and Pa ^ Pti- (31) 

The detunings of central frequencies of CFM and PFM 
are given by A = lu^ — uj^. and 5 = lo^ — lo^, respectively. 
In Eq. ( ^0| ) we have presented the equations for the rel- 
evant density matrix elements only. However in all the 



where npara is the refractive index of the background 
parahydrogen. Therefore the refractive index of the 
medium averaged over the transverse distribution of the 
PFM, using the definition in Eq. (pT|), can be written as 



Sa^pirW^ (J2&{r)]£p{r)dr 



1/2 



(36) 



However assuming that the susceptibilities Xortho and 
Xpara are small, we can rewrite the above as 



''para 



g /f g;(r)[Rea26(r)]gp(r)dr 
2 r|fp(r)Pdr 



(37) 
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Using this rim and the same procedure as described in 
Sec. Ill and IV, we calculate Pp and hence the group 
velocity of the PFM inside the thin fiber surrounded by 
ortho-Il2 doped solid para-Il2 crystal: 



1 



-(/)e)2(l + 2,/)pa) 



2cGl 



--(-/ 



db 



(38) 



assuming the ground state mixing rate 2r = 0. Here 
the parameters Go, ^, 6, (t>p, 4>c are same as defined ear- 
lier. For simplicity we have derived the above equation at 

= and use this condition in all further calculations. 
This analytical result agrees well in the order of magni- 
tude estimation with the numerical result presented in 
the next section. 

In the following section we present the numerical re- 
sults and discussions on various aspects of slow light 
propagation of dressed mode propagating inside the thin 
fiber with orthohydrogen doped in parahydrogen crystal 
as the EIT medium. 



VI. NUMERICAL RESULTS AND DISCUSSION 

The propagation of the probe field inside the thin fiber 
is numerically carried out using the two-dimensional (2D) 
beam propagation method (BPM) [ p2| that solves the 2D 
scalar equation 



+ -g^+n'^{x,z)klAp^O, 



with 
Ap 

£1 



IS 



as the only transverse component, 
the amplitude of the probe field; 



(39) 



Here 



Ap{x;uj)e'''p 



ikjjZ — iiut 



duj-\-c.c. used for initiation of the 



probe field propagation which eventually evolves to the 
fiber mode £p shown in Fig. 4(a). The solution to the 
above equation is obtained by discretizing space and as- 
suming the medium to be a combination of short sec- 
tions of homogeneous medium, having a constant refrac- 
tive index n, separated by lenses at the grid points that 
incorporate the effect of the actual refractive index of the 
medium n{x, z) at the given (x, z). Thus the field propa- 
gation is carried out in two steps: (1) homogeneous step: 
The homogeneous part of the propagation was carried 
out by Fourier decomposition in terms of plane waves; 
(2) lens step: In this step a phase shift is introduced to 
incorporate the effect of n{x,z). We have included fol- 
lowing features into 2D BPM to describe the dynamics of 
our system: (a) a complex refractive index Um {x) is used 
to account for the absorption in the resonant medium; 
(b) the averaged refractive index n for the homogeneous 
step is determined adaptively during the propagation of 
the field using 



^^j^A*p{x)n'^{x)Ap{x)dx 
!%\Ap{x)\Hx ' 



(40) 



thin fiber 




-4.0 -2.0 0.0 2.0 4.0 

X (in nm) 

FIG. 4: Distribution of the (a) PFM (scaled), (b) the real 
part of refractive index, and (c) the imaginary part of the re- 
fractive index that represent the absorption, along the trans- 
verse direction. 
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FIG. 5: Variation of slope of the real part of the refractive 
index with respect to the PFM frequency in the transverse 
direction. 



and (c) the energy conservation \Ap{x)\'^dx]/ dz = 

is used to compensate the numerical losses in the ho- 
mogeneous steps. 

The numerical evaluation of the dressed fiber mode 
PFM, shown in Fig. ^(a) involves the following steps. 
First step is: to determine the transverse distribution 
of the CFM £c{x) using Eqs. (^) and ( ^ with n„i = 
n-paxa = 1.12. We use this £c{x) in Eq. ( |33| ) and in all 
further calculations, and evaluate the response of the 
medium to the PFM in the presence of the CFM. We 
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FIG. 6: The dispersion (solid line) and the absorption (dot- 
ted line) characteristics of the dressed fiber mode when the 
thin fiber is surrounded by ortho-H2 doped para-H2 crystal 
as the EIT medium. Here all the parameters are same as in 
Fig. and the PFM detuning 5 is scaled with 7. 

initiate the propagation of probe field with a Gaussian 
field having FWHM 2a and numerically propagate us- 
ing the 2D BPM as described above. In few steps, the 
Gaussian probe evolves to the PFM E^ix) as shown in 
Fig. 0(a). The distribution of refractive index and ab- 
sorption along the transverse direction are shown in Fig. 
^(b) and (c) respectively. Clearly the behavior of the 
PFM is mainly determined by the region close to this 
fiber (for < 2 /im), as S-(x) is very small outside this 
region. The absorption is large beyond this region as the 
control field reduces and hence EIT condition does not 
hold beyond this region. The parameters used here are: 
iV « 1.3 X 10^7 m-3, deff ~ 7.3 x IQ-^^ Coul.m, 27 = 30 
kHz, r = 1.17 X 10~^7, the resonant Aq = 2.4 /im, 
(5 = — O.OOI7 and {£c)max = ^'^/dcB, the correspond- 
ing laser intensity is « 0.6 W/cm^. The above value of 
iV corresponds to 5% doping of ortho-Il2 molecule in a 
matrix of para-H2 crystal. However since the inhomo- 
geneous width in a solid is much larger compared to its 
natural linewidth, to achieve a good EIT condition, a 
larger intensity is required to overcome the inhomoge- 
neous width Using an effective width of 27 = 20.03 
MHz - where the additional 20 MHz is to account for a 
possible inhomogeneous broadening (27inh) in the crystal 
- the intensity is estimated to be ss279 kW/cm^. This in- 
tensity for our configuration corresponds to a laser power 
of « 19.7 mW only (assuming the beam diameter of the 
PFM of the thin fiber to be ~3 /zm). Further some early 
experimental observations have predicted possibility of a 
much smaller inhomogeneous width for a hydrogen crys- 
tal prepared under very good experimental condition - 
that would reduce the power requrcment for the control 
field by square times the the amount by which inhomo- 
geneous width is reduced; e.g. if inhomogeneous width is 
20 times smaller, the power required will be reduced by 
400 times. Here for the approximate estimation, 7inh is 
assumed to be Lorentzian and added to the natural 
linewidth to get the effective 7. In Fig. we show the 



variation of the slope of the refractive index along the 
transverse direction. It is clear from Fig. I that the 
slope is positive only in the region close to the fiber, 
and hence favorable for the slow group velocity. For 
\x\ > 2 /im the slope decreases and becomes negative 
for |x| > 4 /im. However, the contribution from the later 
part is extremely small as the PFM £p{x) has very small 
value beyond > 2 /im. Further £p{x) is used to deter- 
mine the complex rim using Eq. ( p6|) a nd hence Bp could 
be determined using Eqs. (0) and (|12|). In Fig. o we plot 
the scaled Pp and Im n„i as function of the PFM detun- 
ing S. The corresponding group velocity Vg of the PFM 
«44.1 m/sec. The resulting group delay of the PFM is 
as large as ^1.13 /xsec for a medium length of 50 /im 
(along the direction of propagation). For completeness 
we note that the group velocity of a weak probe in free 
space, passing through the bulk solid sample of ortho- 
hydrogen doped in a parahydrogen matrix, is estimated 
to be Vg ~ 52.95 m/sec. We have used a control field 
strength equal to £c{r = a); all other field and medium 
parameters are taken to be the same as in the above. It 
may be noted that the group velocity reduction in a fiber 
is larger than that in a bulk EIT medium (without fiber), 
which is due to following reason: group velocity of probe 
pulse is directly proportional to the control field inten- 
sity. In case of a fiber, the control field itensity decreases 
along the transverse direction causing the tail part of the 
mode to move extremely slowly, which reduces the group 
velocity of the fiber mode drastically. However in the 
case of a bulk sample, the intensity is assumed to be con- 
stant along the transverse direction. Therefore with the 
same peak value for control field, the probe field pulse 
velocity is slower in the the fiber compared to that in the 
bulk medium. 

It may further be noted that f3p is a constant of prop- 
agation only if we assume that the absorption of both 
pump and probe are negligible. Because a change in 
pump or probe intensity distribution could change rim 
and hence f3p. Validity of this approximation lies in the 
fact that we are working under EIT condition Im hm — 
around S = - which, however, is good only in the re- 
gion close to the fiber. Moreover the probe field ampli- 
tude drops rapidly outside this region. Hence the probe 
absorption and dispersion characteristics are primarily 
determined by the region close to the fiber, where EIT 
condition is satisfied. 

Apart from the usefulness of our configuration to re- 
duce the group velocity of the fiber mode, it offers fol- 
lowing advantages: a fiber mode, of the original (undis- 
torted) fiber with large width at the input - when prop- 
agates through the thin (tapered) region, the width is 
considerably reduced. Thus the tapered fiber acts fike a 
focusing element and hence a large part the field energy is 
available in the evanescent-like part of the thin fiber mode 
- that could be useful to access weaker atomic/molecular 
transitions such as the Qi(l) transition of ortho-H2. 
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VII. CONCLUSION 

In summary, we have demonstrated that the group ve- 
locity of a fiber mode propagating inside a thin single- 
mode optical fiber can be reduced drastically utilizing the 
EIT. We have shown that due to the strong interaction of 
the evanescent-like part of the fiber mode with the EIT 
medium surrounding the fiber, the fiber mode picks up 
the dispersion and absorption properties of the medium. 
We have analyzed in detail how a dressed fiber mode 
can be prepared and the slow fiber mode propagation 
can be realized when the fiber is surrounded by an EIT 
medium, both numerically and with approximate analyt- 
ical results. We have introduced orthohydrogen doped in 



parahydrogen crystal as a new solid EIT medium and 
demonstrated in detail that using this as the medium 
surrounding the fiber, we can reduce the group veloc- 
ity of the fiber mode to ~44 m/sec for the chosen set of 
CFM parameters. From the analytical result we also pre- 
dict the possibility of stopping the PFM by dynamically 
switching off the CFM. 
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